An ARPES Study of the Electronic Structure of the Quantum Spin Liquid EtMe3Sb[Pd(dmit)2]2 
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The electronic structure of a quantum spin liquid compound, EtMe3Sb[Pd(dmit)2]2, has been studied with 
angle-resolved photoemission spectroscopy, together with two other Pd(dmit)2 salts in the valence bond solid 
or antiferromagnetic state. We have resolved several bands that have negligible dispersions and fit well to the 
calculated energy levels of an isolated [Pd(dmit)2]2 dimer. EtMe3Sb[Pd(dmit)2]2 being a Mott insulator, its lower 
Hubbard band is identified, and there is a small gap of ~50 meV between this band and the chemical potential. 
Moreover, the spectral features exhibit polaronic behavior with anomalously broad linewidth. Compared with 
existing theories, our results suggest that strong electron-boson interactions, together with smaller hopping and 
on-site Coulomb interaction terms have to be considered for a realistic modeling of the organic quantum spin 
liquid systems like the Pd(dmit)2 salt. 
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PACS numbers: 71.36.+C, 79.60.Fr, 71.20.Rv 

Quantum spin liquid (QSL) is an exotic state of quantum 
matter, which has profound implications for magnetism and 
high temperature superconductivity 1 1 ]. For two-dimensional 
(2D) spin- 1/2 systems in the presence of antiferromagnetic 
(AFM) interactions, QSL might be realized if the geometrical 
frustration is so strong that it completely destroys the long- 
range magnetic order even at zero temperature, giving rise to 
a large degeneracy of the ground states. Although the QSL 
state on a 2D triangular lattice was first proposed in 1973 |i2|, 
it is hard to achieve experimentally, because the finite three- 
dimensional (3D) magnetic interactions often help drive the 
system into certain ordered states. Recently, some promising 
QSL candidates have been reported 11314 1211. for example, k- 
(BEDT-TTF)2Cu2(CN)3 f?) and EtMe3Sb[Pd(dmit)2]2 (A on 
2D triangular lattices, BaCu3V208(OH)2 |7] on a Kagome lat- 
tice, Na4lr308 181] on a hyper-Kagome lattice, and Ho2Ti207 
on a pyrochlore lattice [9] . They all show no sign of long- 
range order down to very low temperatures. 



EtMe3 Sb [Pd(dmit)2] 2 (dmit = 1 , 3 -dithiole-2-thione-4,5 - 
dithiolate, Me=CH3, Et=C2H5), as shown in Fig. [TJa), 
provides a physical realization of a frustrated spin- 1/2 QSL 
system on a triangular lattice ||4l [isl uM- In fact, a series 
of layered organic salts X[Pd(dmit)2]2 are Mott insulators 
at ambient pressure OlSll . where X is a nonmagnetic mono- 
valent cation such as Et2Me2P, EtMe3Sb, EtMe3P, etc. The 
Pd(dmit)2 units are strongly dimerized, forming a spin- 1/2 
unit [Pd(dmit)2]2 [Fig. db)]. The recent first-principles 
calculations on a series Pd(dmit)2 salts have shown that the 
choice of X would affect the hopping integrals (f, f), and the 
strength of the frustration {J' I J) JmllTll . where the exchange 
interaction terms J' - Jr and J - Jg - Js [Fig. [lib)]. On 
the other hand, the effective on-site Coulomb interaction (U) 
is weakly dependent on the choice of X 11511 . The different 
parameters would then give different ground states lll6lll7ll . 
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FIG. 1: (color online) The structure of EtMe3Sb[Pd(dmit)2]2. (a) 
A view from the b axis. The crystal structure is monoclinic, and 
consists of two-dimensional Pd(dmit)2 layers separated by layers of 
cation EtMcsSb. (b) An in-plane sketch of the spin arrangement on 
the Pd(dmit)2 layer. Each lattice (encircled by solid lines) stands for 
strongly dimerized [Pd(dmit)2]2. Three nearest-neighbor exchange 
interactions, Jb, Js and /,-, are similar and estimated to be 220 - 
250 K. (c) The Laue picture of the EtMe3Sb[Pd(dmit)2]2 single crys- 
tal. Because the dimer units in panel a stack along the a+fr direction 
in one layer and along the a-b direction in the next layer |15], the 
crystal shows twofold symmetry. The inset shows the photo of a 
single crystal. 



For example, besides the QSL, Et2Me2P[Pd(dmit)2]2 ex- 
hibits AFM order below ~ 20 K lUl, EtMe3P[Pd(dmit)2]2 



is a valence bond solid or VBS, where two neighboring 
sites form a spin singlet without magnetic order [, 19,1 . and 
Et2Me2Sb[Pd(dmit)2]2 is in a charge ordered state 112011 . 

All the QSL, VBS, and AFM Pd(dmit)2 salts are Mott insu- 
lators, so they are often modelled using the Hubbard model. 
The band calculations give metallic ground states with siz- 
able f IM Int |2ll |22|]. Then Uke for the cuprate, U, or the 
ratio of UJt, is critical in determining different ground states 
here 1231 12411 . It has been argued that by reducing U jt, the 
system goes from a phase with long-range spiral magnetic or- 
der, to a spin liquid state, and then to a metallic phase ll23ll . 
That is, the spin liquid state on a triangular lattice is stable 
over a narrow range on the verge of metal-insulator transition, 
suggesting that the electron itinerancy may play a role in re- 
alizing spin liquid state 123142511 . More specifically, for the 
EtMe3Sb[Pd(dmi t)2l2 Q SL compound, its t was estimated to 
be ~ 30-55 meV ila,ll7ll22ll . and U has been estimated to be 
about 0.7 eV 12211 . However, experimentally, little is known on 
the basic electronic structure of a QSL, such as its band struc- 
ture and Mott gap (or U), and whether the Hubbard model is 
sufficient. 

In this paper, we report the detailed electronic structure of 
EtMe3Sb[Pd(dmit)2]2 measured with angle-resolved photoe- 
mission spectroscopy (ARPES). We find three broad spec- 
tral features with negligible dispersion within the first L5 eV 
binding energy range, which match well with the calculated 
energy levels of the isolated dimer The weak dispersion and 
broad lineshape resemble the electronic structure in polaronic 
systems, indicating strong interactions between electrons and 
certain bosonic modes in the system. The feature near the 
chemical potential (yu) is recognized as the lower Hubbard 
band (LHB), and the Mott insulator character is evidenced 
by the lack of spectral weight within ~50 meV of the chemi- 
cal potential. QSL compound is found to be more insulating 
than VBS and AFM Pd(dmit)2 salts. Our results indicate that 
the experimental t and U terms of EtMe3Sb[Pd(dmit)2]2 are 
much smaller than what were previously considered, which 
naturally explain its various anomalous properties. We sug- 
gest that the Hubbard model is not sufficient for describing 
EtMe3Sb[Pd(dmit)2]2 or the organic QSL systems in gen- 
eral, instead, a Holstein-Hubbard model with strong electron- 
phonon/magnon interactions is more appropriate. 

The ARPES measurements were conducted at Beam- 
line 7 of Ultraviolet Synchrotron Orbital Radiation Facility 
(UVSOR). The overall energy resolution was set to 20 meV 
or better and the typical angular resolution is 0.5°. Fine single 
crystals of EtMe3Sb[Pd(dmit)2]2, EtMe3P[Pd(dmit)2]2, and 
Et2Me2P[Pd(dmit)2]2 were synthesized by an aerial oxidation 
method llI5[l . pre-oriented by Laue diffraction, then cleaved 
in situ and measured under ultrahigh vacuum of 6 x 10"" 
mbar During measurements the temperature was kept at 70 K 
to avoid the onset of charging at 40 K. Most importantly, we 
employed low energy photons as the excitation light source 
and kept the samples exposed under low photon flux to reduce 
the radiation damage, so that the aging effects are negligible 
in the data. 
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FIG. 2: (color online) Photoemission data of EtMe3Sb[Pd(dmit)2]2. 

(a) The two-dimensional Brillouin zone of EtMe3Sb[Pd(dmit)2]2. 

(b) The normal-emission photoemission spectra of 
EtMe3Sb[Pd(dmit)2]2 excited by p, s, left-handed circular (C-I-) 
and right-handed circular (C-) polarized 8 eV photons. The inset 
shows the normal emission spectra collected with p polarized 8, 
15, 21 eV photons, (c) The photoemission intensity along the (0, 
-;r/a)-(0, ;r/a) line as indicated in panel a. Here in order to highlight 
the main spectral features, the original data has been deducted off 
the Shirley background obtained in Fig.|3la) to exclude the influence 
of secondary electrons, (d) ARPES spectra of EtMe3Sb[Pd(dmit)2]2 
acquired at momenta #1 - #15 as coded in panel a. The photoe- 
mission data in panels c and d were taken with p polarized 8 eV 
photons. 



The single crystal of EtMe3Sb[Pd(dmit)2]2 is usually a 
hexagonal piece with typical size of 0.6 x 0.6 x 0.05 mm^ as 
shown in the inset of Fig. [lie). The unit cell is monoclinic as 
shown in Fig.lTJa) and the space group is C2/c. More detailed 
description of the crystal structure could be found in Ref. llI5ll . 
The high sample quality was confirmed by the typical Laue 
pattern [Fig. [He)], which manifests the twofold symmetry of 
the crystal structure. 

The normal-emission spectra from the F-point of 
EtMe3Sb[Pd(dmit)2]2 BriUouin zone [Fig. I^a)] are dis- 
played in Fig. |2|b). In order to make sure all states are 
detected, various photon polarizations were exploited, includ- 
ing the in-plane (/?), perpendicular-to-plane (s), left-handed 
circular and right-handed circular polarizations. Here, the 
"plane" is defined by the incident photon beam direction 
and the outgoing photoelectron trajectory, so that s and p 
polarizations are sensitive to only odd and even electron 
wavefunctions with respect to this mirror plane jlO, U.l\. We 
found three broad features within the 1.5 eV binding energy 
range, and they are most pronounced in the p-polarization 
data. Moreover, their binding energies are robust against the 
variation of the photon energy as illustrated in the inset of 
Fig. |2|b), which is a manifestation of the 2D character of the 
system. Hereafter, we only show the data taken with p polar- 



ized 8 eV photons, which give the most pronounced spectral 
features. Fig. |2c) illustrates the photoemission intensity 
along (0, -;7r/a)-(0, ;7r/a) after deducting a Shirley background 
as show in Fig. |3ja), which shows spectral features with 
negligible dispersion. Furthermore, as displayed in Fig. |2ld), 
all the photoemission spectra taken at various momenta in the 
first Brillouin zone show three specific spectral features as 
marked by the dashed lines. However, no obvious momentum 
dependence of the peak positions is observed, indicating the 
states are considerably localized in EtMe3Sb[Pd(dmit)2]2, 
contradicting the pronounced dispersions up to 400 meV 
predicted along (0, -7r/a )-(0, 7r/a) and (-7r/b, 0)-(;r/b, 0) in the 
calculations I16l ll7ll22ll . 

The photoemission data of EtMe3Sb[Pd(dmit)2]2 show no 
density of state (DOS) at the chemical potential [Fig. Oa)], 
as expected for a Mott insulator Recent torque magnetom- 
etry measurement on EtMe3Sb[Pd(dmit)2]2 revealed an ex- 
tended quantum critical phase, in which low-lying gapless 
spinon excitations behave like the elementary excitations in 
paramagnetic metals with Fermi surface |'28'|. We here do 
not observe any projection of the spinon Fermi surface in 
the charge channel. Furthermore, we found that the pho- 
toemission spectrum is best fitted with three Gaussian com- 
ponents, plus a Shirley-type background in order to account 
for secondary electrons [Fig. EJa)]. To compare with the ex- 
periment, we have calculated the energy levels of the iso- 
lated Pd(dmit)2 monomer and [Pd(dmit)2]2 dimer in a super- 
cell using their structures in EtMe3Sb[Pd(dmit)2]2 with the 
all-electron full-potential linearized augmented plane wave 
(FLAPW) method. The resulting energy levels are almost the 
same as those calculated in Ref. Il29tl for another [Pd(dmit)2]2 
salt. As illustrated in Fig. Ob), the interactions between 
the two Pd(dmit)2 molecules in an isolated dimer will split 
their respective HOMO (highest occupied molecular orbital) 
and LUMO (lowest unoccupied molecular orbital) states into 
bonding and antibonding states, which are labeled as Hb, Ha, 
Lb and La self-explanatorily. Based on the stoichiometry, 
there are five electrons, resulting in a half-filled Ha band. It is 
also the fully filled LHB of this Mott insulator Based on the 
fitting, the three Gaussian components located at 0.31 + 0.02, 
0.67 ± 0.02, 1.20 + 0.02 eV and labeled as peakl, peak2, and 
peak3, which correspond to the //a/LHB, Lb, and Hb bands, 
respectively [Fig. lUb)]. The agreement between the experi- 
ment and calculation is reasonably good, especially consid- 
ering that there are small difference between the intra-orbital 
and inter-orbital onsite Coulomb interactions [22J. Moreover, 
the spectral weight of peak2 and peak3 are comparable, while 
peakl is considerably weaker, which is expected as the LHB 
is only filled with one electron. 

The localized electronic structure is not only manifested 
by the negligible dispersion, but also by the broad linewidth. 
The fitted full width at half-maximum (FWHM) of peakl, 
peak2 and peak3 are 0.23, 0.40, and 0.48 eV, respectively. 
In some calculations, there are a set of bands around 1.2 eV, 
which may cause the possible broad lineshape of peak3, but 
the spectral width of peak2 is much broader than the calcu- 
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FIG. 3: (color online) (a) The normal-emission valence band spectra 
of EtMe3Sb[Pd(dmit)2]2 is fitted with three Gaussians plus a Shirley- 
type background, (b) Left: The centroid of the peaks obtained from 
the fitting. The black arrows represent the electron spins. Right: 
The calculated energy levels of the isolated Pd(dmit)2 monomer 
and [Pd(dmit)2]2 dimer in EtMe3Sb[Pd(dmit)2]2. The antibonding 
LUMO, antibonding HOMO, bonding LUMO, and bonding HOMO 
energy levels are labeled as La, Ha, Lb, and Hb, respectively, (c) 
A schematic spectral function of a single electron coupled to a bath 
of bosons of frequency cog in the weak coupling and strong coupling 
conditions, (d) The spectrum near the chemical potential is enlarged, 
and a gap of 50 meV is hereby defined. 



lated bandwidth of 0.1 eV il6lll7ll22ll . The broad lineshape 
that could be fitted with the incoherent Gaussian peaks re- 
sembles those of polaronic systems, where they are caused 
by strong electron-phonon coupling [30-32]. As shown in 
Fig.[3jc), the quasiparticle weight is usually vanishingly small 
at the low-energy edge of the spectrum, and the broad spec- 
trum is an envelope of many individual peaks due to strong 
electron-boson interactions 1.33,1 . For a flexible organic com- 
pound with strong spin fluctuations, strong electron-phonon 
interactions and/or electron-magnon interactions are expected 
1 3411 . Therefore, the charge carriers in QSL could be regarded 
as polarons, which are heavily dressed by phonons (lattice de- 
formations) or magnons. In addition, it has been shown that 
the centroid of the incoherent spectral distribution would fol- 
low the bare band dispersion in the absence of strong corre- 
lations [35, 36]. Therefore, our data indicate that t is at most 
a few meV, much smaller than that predicted in the previous 
calculations il6lll7ll22ll . It shows that dimers are rather iso- 
lated in the crystal, which is consistent with the fact that the 
centroid positions of the spectral features match the energy 
levels of the isolated dimer in Fig. |3lb). 

Taking the lower edge of the LHB band as the coherent 
quasiparticle location of the polaronic lineshape [Fig.|3jc)], a 
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FIG. 4: (color online) (a) The normal-emission spectra of 
Et2Me2P[Pd(dmit)2]2 (AFM compound) and EtMe3P[Pd(dmit)2]2 
(VBS compound) are fitted with 3 Gaussian peaks plus a Shirley- 
type background, respectively. Data were taken with p polarized 
8 eV photons, (b) The F-point spectra of the three Pd(dmit)2 salts 
in the vicinity of the chemical potential, (c) The schematic diagram 
for the experimental (exp.) and calculated (cal.) energy-levels of 
these materials. The calculated Ha bands are individually aligned to 
the experimental positions of their peakl's. The unit is eV. 



gap could be estimated as in Fig.[3td). It is about 50 meV, or 
smaller, because the DOS gradually decreases to zero at the 
chemical potential. If QSL appears only in a narrow region of 
the U/t axis 023112411 . the small hopping integral t here requires 
the actual U to be much smaller than the previously predicted 
0.7 eV for EtMe3Sb[Pd(dmit)2]2 lH. In fact, this is naturally 
expected when strong electron-phonon interaction is consid- 
ered, polarization of the lattice would significantly reduce the 
effective on-site Coulomb repulsions. Assuming U ~ S t ll23ll . 
U is probably in the order of 50 meV, which explains the 
rather accurate calculated energy levels, since the small dif- 
ferences between inter-orbital and intra-orbital Coulomb in- 
teractions will not give much deviation. Consequently, the 
upper Hubbard band should not be far away from the LHB in 
EtMe3Sb[Pd(dmit)2]2, which naturally explains the low exci- 
tation energy of about 40 meV estimated from the resistivity 
1 37t] and the low charging onset temperature of 40 K in our 
photoemission study. In other words, our results do suggest 
that EtMe3Sb[Pd(dmit)2]2 is on the verge of a metal-insulator 
transition, in agreement with the theory [23. .24.1 . Consistently 
for /c-(ET)2Cu2(CN)3, which is another organic QSL on a tri- 
angular lattice, the charge gap is very small in optical con- 
ductivity measurements (the exact number is hard to extract, 
probably in the 15 ~ 100 meV range) 113 811 . indicating the elec- 
tronic structure characters found here might be fairly generic. 
For comparison, the electronic stiTicture of 
Et2Me2P[Pd(dmit)2]2 (AFM) and EtMe3P[Pd(dmit)2]2 
(VBS) are investigated, whose ARPES spectra with neg- 
ligible dispersion are qualitatively similar to those of 



EtMe3Sb[Pd(dmit)2]2. The typical spectra for these com- 
pounds in Fig. Ua) demonstrate the insulating behavior of 
these Pd(dmit)2 salts, while the ARPES spectral features in 
AFM and VBS materials are broader than those of QSL, 
and their energy separations are displayed in Fig. EIc). The 
calculated energy levels of these Pd(dmit)2 salts are also 
presented in Fig.HJc), which show reasonable agreement with 
the experiments. Near the chemical potential, the spectral 
weight of the QSL is further suppressed from the chemical 
potential than the other two, as clearly illustrated in Fig. lUb). 
For the QSL compound, peakl is located at 0.31 + 0.02 eV 
with narrower width than those of VBS and AFM compounds, 
making the DOS of EtMe3Sb[Pd(dmit)2]2 more suppressed 
near the chemical potential. This explains the relatively lower 
excitation energies of AFM (~ 26 meV) and VBS (~ 27 
meV) compounds than that of QSL compound, as estimated 
from their resistivity data lHH Bll IH. Our ARPES data 
could not provide information on the evolution of t, t' , and 
U in these compounds, because of their small amplitudes. 
However, we do observe certain vaiiations in the energy level 
positions and spectral widths, which suggest the variations in 
the dimer structure and electron-boson interaction strength, 
respectively. Whether and how these tip the balance between 
different ground states require further investigation. 

To summarize, we have presented the experimental elec- 
tronic structure of EtMe3Sb[Pd(dmit)2]2, a QSL, for the first 
time, which naturally explains many phenomena observed 
in this compound. Moreover, we found that Pd(dmit)2 salts 
in the QSL, VBS, and AFM states are Mott insulators with 
small gap and vanishing spectral weight at the chemical po- 
tential. Contrary to previous theoretical picture, we found 
that these systems are characteiized by negligible hopping and 
broad polaronic lineshapes. It indicates that strong electron- 
phonon/magnon interactions have to be included when de- 
scribing these systems. Our results suggest that a Holstein- 
Hubbard model with strong electron-boson coupling, instead 
of the Hubbard model, is more appropriate for a realistic mod- 
eling of the organic QSL's and related compounds. 
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